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The Polyakov loop variable serves as an order parameter to characterize the confined and de- 
confined phases of Yang-Mills theory. By integrating out the vector fields in the SU(2) Yang-Mills 
partition function in one-loop approximation, an effective action is obtained for the Polyakov loop 
to second order in a derivative expansion. The resulting effective potential for the Polyakov loop is 
■ capable of describing a second-order deconfinement transition as a function of temperature. 
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^ ^ ^ Introduction. Confinement of quarks and gluons into color-singlet clusters is the most important characteristic of 
strong interaction physics. The phenomenon has been observed in calculations based on the lattice formulation of 
^/^ Yang-Mills theory, and diverse physical mechanisms for the emergence of a confining potential between color sources 
have been proposed, among others, the dual Meissner effect picture of 't Hooft and Mandelstam random fluxes 
Q , the stochastic vacuum jj] , the leading-log model of Adler |^ , and dual QCD |^ . 

On the other hand, if one generalizes to finite temperature, it is believed (and borne out e.g. by lattice calcula- 
tions that there exists a transition to a phase in which color sources become deconfined. An order parameter 
I> ' distinguishing the two phases in pure gauge theory is the Polyakov loop , Q 
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L{x) = -rj— tr Texp I i I dxQAQ{xQ,x) , (1) 




■ where /? is the extension in time direction of the (Euclidean) space-time manifold under consideration. It can be 
identified with the inverse temperature of the ensemble in which the expectation value of L is taken. Aq denotes 
, the time component of the gauge field, T time ordering, and Nq the number of colors. The negative logarithm of 
' the expectation value of the Polyakov loop can be interpreted |^ as the free energy associated with a single static 
, color source in the fundamental representation of the gauge group. Vanishing of (L) indicates infinite free energy, 
i.e. confinement; finite {L) indicates deconfinement. There is a loophole to the former interpretation: The free 
(-H I energy may be infinite due to reasons not connected with the infrared phenomenon of confinement, e.g. an infinite 
self-energy, which is an ultraviolet phenomenon. One should thus take care to separate infrared and ultraviolet effects 
, in the evaluation of (L). The behavior of (L) is intimately connected to a discrete symmetry of Yang-Mills theory 
. called center symmetry . This symmetry leaves the Yang-Mills action invariant but changes L by a phase q'^^^/^c _ 
$^ ' Therefore, if the symmetry is realized, one automatically has (L) = 0, i.e. confinement; if the symmetry is broken, one 
5^ , has (L) ^ 0, deconfinement. Which scenario is realized can be read off by calculating the effective action governing the 
Polyakov loop order parameter and interpreting it in the spirit of Ginzburg-Landau theory Such considerations 
have been quite successful in the lattice formulation ]To|- ]T^ , where, however, as always the continuum limit represents 
a problem. A continuum calculation was carried out in | |13[ , which was only able to capture the deconfined phase (see 
however also JT^ , in which a model built on ideas familiar from the Savvidy vacuum | |l5|] is suggested) . 

In this letter, a continuum evaluation of the effective action is exhibited for the case of SU(2) color which goes 
significantly beyond previous work and in contradistinction to describes both the confined and deconfined phases. 
Specifically, centrifugal barriers in kinetic terms which were not considered in substantially modify the effective 
potential for the Polyakov loop order parameter. The success of this calculation gives rise to the hope that the confined 
phase can be understood essentially from the metric properties of the gauge group without extensive recourse to details 
of the dynamics. 

Approximation of the SU(2) effective action. A particularly convenient gauge for discussing the behavior of 
the Polyakov loop variable (0) in the case of SU(2) color is @- 

9oAo = and = diag (ao/2, -ao/2) . (2) 

In this gauge, the Polyakov loop takes the simple form L = cos/3ao/2 and center symmetry relates oq and 27r//3 — ao 
if ao S [0, 2tt/P]. Thus, in order to gain insight into the dynamics of the Polyakov loop, one should integrate out in 
some approximative scheme the vector fields Af in the SU(2) Yang-Mills partition function [H 
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Z = J [Dao] detfiiao) J [i?4°]e-^--['^'"^?l (3) 

(without much loss of generality, only the case of vanishing vacuum angle Q is considered here). The measure det /i(ao) 
denotes the Cartan part of the SU(2) Haar measure for uq at every space point. 

The approximation to be used here is the following: Terms in Sym of higher order than quadratic in the vector fields 
are dropped, and the resulting functional determinant is calculated in a gradient expansion up to terms quadratic in 
derivatives of gq. Previous work only considered constant oq. It will become apparent that centrifugal barriers in 
the kinetic terms lead to a substantial modification of the effective potential. The effective action for oq in the above 
approximation can be written as 

Seff[ao] = / (fixdiaodiaQ — Indet^(ao) — lndet~^/^A , (4) 

where the last term stems from the Gaussian integration over the vector fields A^. The operator A becomes block- 
diagonal in color space if one changes to the variables {Aj ± iAf)/\/2, A^. Then one has two blocks 

^fj = -^^3^^ + + % i-tdo ± ao)^ (5) 

and a third block independent of gq which gives only an irrelevant constant. 

Evaluation of the determinant. Concentrating on the operator A~, its functional determinant can be expressed 
in the form 1191 



lndet"^/^A" = - tttttt / dtf^'^f] / d^xtYK(xo,x;xo,x;t) 
2 as r{s) Jo J 



(6) 

s=0 



with the heat kernel 

Kij{xo,x]yo,y]t) ^ (^e^*^ _ S^{x - y)d{xQ - yo) , (7) 



where A is evaluated at, and acts on, the coordinates x and xq. In the limit (x — y) — > 0, one can expand 



ao{x) = ao{y) + {x - y)idiao{y) + -{x - y)i{x - y)jdid.jao{y) . (8) 



Higher orders in derivatives of ao are neglected, as mentioned above. Upon inserting Fourier representations for the 
(5- functions in ^ and commuting exp(— tA~) through to the right, one obtains the representation 

1 oo „ 

K,^{xo,x;yo,y;t) = -^ ^ / ^3fcg^(;c(x-y)+2.K-o-yo)//3)e-tfe=g-t(2W^-ao(v))^ (^^-(a+s)^^ (9) 



with 



B^j = t{~6ij{2ikidi + d^) + ihdj + ik^di + didj + - y)mCm + Sij{x - y)m{.x - y)nDmn) (H) 
C™ = 2(ao(2/) - 2TTr/l3)drnao{y) (12) 
Dmn = {ao{y) - 2Trr/P)dmd„ao{y) + dmaQ{y)d„aQiy) . (13) 

A and B are to be interpreted as operators acting to the right on unity. Now one can expand in powers of B, since 
only terms up to quadratic in derivatives of oq are to be retained. One needs terms up to fourth order in B. Since A 
and B do not commute, the expansion 



f dsie-(i-'*i)^Be-''i^ (14) 

+ [ ds2 [ dsiSie-''^-''^'^Be-''^^-'^^'^Be-'''^'^ - ... 
Jo Jo 

must be used. In the limit a; — y — *■ 0, only terms survive in which all factors (a: — y)i are removed by derivative 
operators. The matrix algebra is carried out using {exp(—wA))ij = Sij + (e'xjp{wtk'^) — l)kikj /k"^ . Furthermore 
evaluating the integrals over the Si in (fil), one arrives at 
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If °° 
/3(27r)3 J 



2e-^^' + 1 



+ (2i(l + e-**^ )lk'^ - 4(1 - e"**^ )lk^) ■ {h{D,j - CiCjt/A)kj - k^D,, + k^C^dt/A) 
+e-*''\8hDijkjt^/3~ hdCjkjt* ~2Duf + 2CiC,t^/3)\ . (15) 

In order to regularize the momentum integrations in the uhraviolet, one needs to introduce a cutoff A on |fc|. Sub- 
sequently expanding the terms in the square brackets in powers of t, one can do the i-integration called for in (|^) 
using dt e~*'°'t^~'^^'^ — r(s + n)/a*"'"". After dividing by T{s) and taking the derivative w.r.t. s at s = 0, one 
can resum the expansions (these are derivatives and integrals of geometrical series) . After explicitly inserting Da and 
CiCi (cf. ( p^ ) and (^2|)) in order to make all Matsubara frequencies manifest, one can finally carry out the Matsubara 
summations using 



E 



P sinh/?y 

(16) 



(27rr//3 + x)2 + y2 2y cosh f3y — cos Px 

E o-„/i ^ ^-|cot(/3x/2) , (17) 



2-Kr/P-x 2 

and appropriate derivatives and integrals over x and y thereof. One finally arrives at 

Indet" ' A = 2 • -/3 / d^x— — — ^ In sin^ ttc ^--r / j/^ ln(cosh ?/ — cos 27rc) 

2 J P y Q-K^ TT-^p-i J„ 

-,^O.cd.c ■ - H ^ n d_y sinh, 

2/3 \ 3sin TTc 9 cosh /3A — cos 27rc 3 Jq y cosh y — cos 27rc 

4/3A 1 - cos 27rc cosh /3A /J^ A^ sinh /3A(2 - cos^ 27rc - cosh /3A cos 27rc) \ 
9 (cosh/3A - cos27rc)2 9 (cosh/3A - cos 27rc)3 J 

where one derivative has been shoveled off terms didiUo by means of a partial integration, the rescaled variable 
c = /3ao/27r has been introduced, and the factor 2 in front comes from adding the contribution of the determinant of 
A+ (it turns out that det A~ and det A+ separately are symmetric in ao). The first term in the square bracket in ( |T^ ) 
exactly cancels the term — In det /i(ao) in the effective action (^ stemming from the Haar measure. This cancellation 
was already observed for space-time independent fields in |l^. In that case, after splitting the vector fields into their 
transverse and longitudinal parts, one sees directly that the longitudinal subdeterminant cancels the Haar measure. 
Interpretation and discussion. The calculation sketched above has yielded an effective action of the form 

S,ff[c] = J d'x l^-^Vic) + ^Wic)d.cd.c^ (19) 



for the field c. V{c) is (up to a factor — /3^) the second term in the square brackets in ([l8|), and W{c) is the term in the 
large parentheses in ( p^ ) (up to a factor —1/2), supplemented by a constant from the tree level contribution in (^. It 
should be emphasized that the preferred mean field value of c cannot be read off from V{c) alone. When evaluating 
path integrals over c, only the class of space-time constant configurations will give its dominant contribution at the 
value minimizing V{c). The weight of all other configurations is influenced by W{c), which may strongly penalize 
values preferred by V{c). In order to properly extract an effective potential from this action, one should transform 
the weight function W away from the kinetic term by the local change of variables u{c) = J^^^ dc' ^JW{c') / p. The 

action in the new variables u then acquires the standard form 5*6// [u] — J d^x {diudiU + Vef f{c{u))), where 



VefMu)) - ^V{c{u)) + ^LJ^i^w{c{u)) (20) 



^3 ' v-v-" ■ i27r2/33 

can now be straightforwardly identified as an effective potential. The second term in (|2^) comes from the Jacobian 
of the change of variables in the path integral. Since one is ultimately interested in the preferred value of c, one 
may minimize this potential directly as a function of c. Note that in the /3A — > oo limit, one recovers from V{c) the 
expression given in for the effective potential, V{c) — > (47r^/3)c^(l — c)^. In this same limit, the term in the large 
parentheses in ( p^ entering W{c) is dominated by its first summand. Figure [| displays the effective potential Fg// 
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as a function of c for various large /3A at a suitably small coupling which is the regime one is ultimately interested 
in when taking the cutoff A to infinity. The value c = 1/2 corresponds to vanishing Polyakov loop order parameter. 
One observes a second-order transition between a regime of larger /3A, where the minimum of the effective potential 
is at c = 1/2, and a regime of smaller /3A, where the order parameter develops a nonvanishing expectation value due 
to spontaneous breaking of the center symmetry in the double-well potential. The precise value of /3A at which the 
transition takes place rises as the coupling constant g falls. Some care should be taken in interpreting the role of g. It 
is possible to adjust g as the cutoff A is taken to infinity such as to obtain finite results for observables (more specific 
remarks on this follow further below). However, one should not expect the renormalization to mirror the behavior 
familiar e.g. from perturbation theory; note that the loop expansion carried out for the vector fields cannot be equated 
with the usual perturbative expansion, since the field c entering the Polyakov loop is of order 1/g. Furthermore, the 
renormalization of g in general also depends on the inverse temperature (3. The authors have not explored the 
possibility of establishing a renormalization group type of equation relating renormalized coupling constants g at 
different /3. At this stage, the effective theory can be used as follows: Fit the temperature dependent running coupling 
g such as to reproduce the temperature dependence of a physical quantity, e.g. the string tension in the confining 
phase, always in the physical limit A —f oo. Then other quantities can be predicted. The authors have tested the 
feasibility of this scheme by deriving the string tension from the Polyakov loop correlator. Its negative logarithm, 
— (1//3) ln(L(i?)L(0)), can be interpreted as the free energy associated with a pair of static sources separated by 
the vector R. On the other hand, for small oscillations of the field c around c = 1/2 in the confined phase one has the 
estimate (L(i?)L(0)) = {u{R)u(0))Tr'^ P /W {c = 1/2). Furthermore, in this approximation one can replace the effective 
action for w by / d^x {diudiU + ji'^a'^v?) with /J^cr^ = {d'^V^f f / dc^) / {2W)\c=i/2- Then the propagator {u{R)u{0)) takes 
the simple Yukawa form (u(i?)u(0)) = e~^'^^ /{WttR), so that a can be directly identified with the string tension. In 
order to reproduce a finite string tension as the temperature approaches zero and a finite deconfinement temperature, 
it turns out that one should scale g{/3/(3cr ^ 1) l/(/3crA)^ and g{/3/(3cr oo) ^ l/(/3crA) where (3cr is the inverse 
deconfinement temperature. 

To conclude, a quite drastic approximation on the dynamics of the vector fields has been shown to nevertheless 
reproduce the confining and deconfining phases of Yang-Mills theory. It appears that the confining phase can be 
simply understood as a consequence of the metric properties of the gauge group, encoded in the centrifugal barriers 
entering the kinetic terms, cf. eq. (^^. The detailed dynamics of the vector fields Ai were largely ignored in the 
present approximation. The picture obtained here for the emergence of the confined phase stands in contrast to the 
Z(2) bubble mechanism advocated on the basis of previous work |]l3|. In the low-temperature regime, there is no need 
for a Z(2) domain structure to restore center symmetry. However, Z(2) bubbles may be relevant for the existence of a 
superheated confining phase above what one would usually consider the deconfinement temperature. The confinement 
mechanism in this phase would be akin to a localization phenomenon. A further comment concerns the formulation 
of high-temperature perturbation theory. There, one usually assumes small fluctuations of the field ao around ag = 0. 
Clearly, this is not correct even in the deconfined phase, since the value ao = is always protected by centrifugal 
barriers; one should therefore expand around a nontrivial value for oq as given by the effective potential (pO|). It is 
tempting to speculate |Q that this is the origin of the infrared problems initially reported by Linde |^l|]. Work is in 
progress to calculate perturbative corrections to the effective potential, to extend the treatment to SU(3) color, and 
to determine how other observables, such as the gluon condensate, behave under the present approximation scheme. 
Also a more detailed treatment of the renormalization is under study. 
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FIG. 1. Effective potential, up to an irrelevant constant shift, as a function of c = /3ao/27r. 
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